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We estimated the Galactic model parameters from the Sloan Digital Sky Survey (SDSS) data reduced for two fields in 
the anti-centre direction of the Galaxy, / = 180°, and symmetric relative to the Galactic plane, b = +45° (north field) 
and b — —45° (south field). The large size of each field, 60 deg 2 , and the faint limiting apparent magnitude, go = 22, 

in ■ 

give us the chance to determine reliable parameters for three components, thin and thick discs and halo, in the north 
and south hemispheres of the Galaxy, except the scalelengths for two discs which are adopted from June et al. (2005). 
Metallicities were evaluated by a recent calibration for SDSS, and absolute magnitudes of stars with 4 < M(g) < 8 
were derived as a function of (g — r)o colour and metallicity. A \ 2 method was employed to fit the analytical density 
laws to the observational-based space densities with the addition constraint of producing local densities consistent with 
those derived from Hipparcos. Conspicuous differences could not be detected between the corresponding Galactic model 
parameters for the thin disc of north and south fields, and our results are consistent with the ones in the literature. The 
same case is valid for the halo, especially the axis ratios for two fields are exactly equal, k = 0.45, and close to the one 
of Juric et al. (2005). However, we revealed differences between the scaleheights and local space densities for the thick 
disc of the north and south fields. The metallicity distribution for unevolved G type stars with 5 < M(g) < 6 shows three 
substructures relative to the distance from the Galactic plane: for z* < 3 kpc, the metallicity gradient for two fields is 
d[M/H]/dz ~ — 0.20(±0.02) dex kpc -1 , consistent with the formation scenario of the thin disc. For 5 < z* < 10 kpc, 
d[M/ H]/dz ~ — 0.03(±0.001) dex kpc -1 for two fields, confirming that the metallicity gradient for the halo component 
is close to zero. However, the tendency of the metallicity for stars with 3 < z* < 5 kpc, corresponding to the transition 
region from thick disc to halo is different. For the north field d[M/H]/dz = — 0.36(±0.12) dex kpc -1 , whereas it is 
half of this value for the south field, d[M/H]/dz = -0.18(±0.01) dex kpc" 1 . The origin of this conflict is probably 
due to the different structure of the thick disc in opposite latitudes of the Galaxy. When we combine these substructures, 
however, we find a smooth metallicity gradient for two fields, — 0.30(±0.04) dex kpc -1 . 
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1 Introduction ued to cover large ranges through the most recent work (Ta- 
ble 1). For example, Chen et al. (2001) and Siegel et al. 
The study of Galactic models and parameters has a long ( 20021 * give 6.5-13 and 6-10 per cent, respectively, for the 
history. Bahcall & Soneira ( [1980] ) fitted observations with relative local space density of the thick disc. However, we 
a two-component Galactic model, namely disc and halo, feel that an appropriate procedure would lead us to expect 
while Gilmore & Reid ( fl"983l ) success was to fit their ob- a smaller range and/or a unique value carrying a small er- 
servations with a Galactic model introducing a third com- ror. In deed, in previous papers (Karaali, Bilir & Hamz aoglu 
ponent, i.e. the thick disc. It should be noted that the third 120041 hereafter KBH; Bilir, Karaali & Gilmore [20061 here- 
component was a rediscovery of the "Intermediate Popula- after BKG) we argued that Galactic model parameters must 
tion II" first described in the Vatican Proceedings review of be mass dependent and, since absolute magnitude provides 
O'Connell ( fl"958l ). Due to their importance, Galactic mod- a reasonable proxy for mass, they must vary from one ab- 
els have been attractive for many research centers: they can solute magnitude interval to the next. Therefore, neglect of 
be used as a tool for uncovering the formation and evolution this mass-or absolute magnitude-dependence leads to large 
of the Galaxy. errors in parameter values, and this probably also explains 
Different research groups have been using different meth- their lar 8 e ran 8 es exhibited by the star-count studies listed 
ods to determine parameter values, and these have contin- ' n Table 1 ■ 
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ues in a crucial way, the present project has been designed 
as follows. 

Our data will consist of homogeneous imaging sets of 
SDSS for stars in several fields in different directions of the 
Galaxy. Using these data, we shall apply a parameter esti- 
mation procedure which is based on the most recent calibra- 
tions: 1) absolute magnitudes will be derived as functions 
of both (g — r)o colour and metallicity and 2) metallicities 
will be determined using the new calibration provided by 
Karaali, Bilir & Tuncel (120051 hereafter KBT). 

The ultimate goal of this project is to derive, in a series 
of papers, more reliable values of Galactic model parame- 
ters, including their dependence on Galactic longitude and 
latitude. 

In Section 2 we describe the SDSS data and their reduc- 
tions, and in Section 3 we present the density laws that we 
shall adopt in the present study. Sections 4 and 5 provide 
the metallicity calibration, the absolute magnitude determi- 
nation, and the evaluation of the density functions for in- 
dividual population components. In Section 6 we estimate 
the stellar metallicities and final conclusions are drawn in 
Section 7. 

2 The Sloan Digital Sky Survey 

The SDSS is a large, international collaboration project set 
up to survey 10 000 square-degrees of sky in five optical 
passbands and to obtain spectra of one million galaxies, 
100 000 quasars, and tens of thousands of Galactic stars. 
The data are being taken with a dedicated 2.5m telescope 
located at Apache Point Observatory (APO), New Mexico. 
The telescope has two instruments: a CCD camera with 30 
2048 x 2048 CCDs in the focal plane and two 320 fiber dou- 
ble spectrographs. The imaging data are tied to a network of 
brighter astrometric standards (which would be saturated in 
the main imaging data) through a set of 22 smaller CCDs in 
the focal plane of the imaging camera. An 0.5m telescope at 
APO will be used to tie the imaging data to brighter photo- 
metric standards. 

The SDSS obtains images almost simultaneously in five 
broad bands (u, g, r, i and z) centered at 3540, 4760, 6280, 
7690 and 9250 A, respectively (Fukugita et al l 19961 Gunn 
et al. [19981 Smith et al . 120021 Hogg et al. 120021 . The imag- 
ing data are automatically processed through a series of soft- 
ware pipelines which find and measure objects and pro- 
vide photometric and astrometric calibrations to produce a 
catalogue of objects with calibrated magnitudes, positions 
and structure information. The photometric pipeline (Lup- 
ton et al. 2001) detects the objects, matches the data from 
the five filters, and measures instrumental fluxes, positions 
and shape parameters. The last allow the classification of 
objects as "point source" (compatible with the point spread 
function), or "extended". The psf magnitudes are currently 
accurate to about 2 per cent in g, r and i and 3-5 per cent in 
u and z for bright (< 20 m ) point sources. The SDSS in Data 
Release 4 (DR4) is almost 85 per cent complete for point 



sources to (u, g, r, i, z) = (22, 22.2, 22.2, 21.3, 20.5), and 
the full-width at half-maximum of the psf is about 1.5 arc- 
sec (Abazajian et al. 2004). The data are saturated at about 
14 m in g, r and i and about 12 m in u and z. Astrometric 
calibration is carried out using a second set of less sensitive 
CCDs in the camera, which allows the transfer of astromet- 
ric catalogue positions to the fainter objects detected by the 
main camera. Absolute positions are accurate to better than 
0.1 arcsec rms in each coordinate (Pier et al. 2003). 

2.1 Data and reductions 

The data are taken from SDSS (DR4) on the WEr£] for two 
fields with the same Galactic longitude (I — 180°) and sym- 
metric relative to the Galactic plane, b — +45° (hereafter 
north) and b = —45° (hereafter south). The areas of the 
fields are equal and rather large, 60 deg 2 , avoiding any fluc- 
tuation in the values of Galactic model parameters. SDSS 
magnitudes u, g, r, i and z are available for 516 279 point 
sources in two fields down to the limiting magnitude of 
g a = 24. The total absorptions, A u , A g , A r , Ai and A z , 
for each point source are taken from SDSS query server, and 
the apparent magnitudes u, g, r, i and z are de-reddened by 
equations in eq. (1) 

u a = u - A u , g = g - A g , r = r - A r 

i = i - Ai,z = z - A z , (1) 

All the colours and magnitudes mentioned hereafter will be 
de-reddened ones. 

Star/galaxy separation was performed by utilizing the 
command prob PSFmag in DR4 WEB page, to provide each 
objects probability of being a star in each filter (parameter 
values of or 1 indicating a galaxy or star, respectively). 
The quality of this separation therefore strongly depends on 
the seeing and sky brightness. In addition to the work cited 
above, we adopted the procedure of Juric et al. {2005 ) in or- 
der to remove the scattered point sources in the (r — i)o — 
(g — r)o diagram. Thus, sources at distances larger than 
0.3 mag from the locus were excluded from the statistics 
(Fig. 1). Such large deviations are inconsistent with mea- 
surement errors, and in most cases indicate a source that 
is not a main sequence star. This requirement effectively re- 
moves hot white dwarfs (Kleinman et al. 2004), low-redshift 
quasars (z < 2.2, Richards et al. 120011) . and white/red dwarf 
unresolved binaries (Smolcic et al. 2004). 

The apparent magnitude histograms for point sources 
and for the star sample are shown in Fig. 2, where the limit- 
ing apparent magnitude for the survey stars can be read off 
as go = 22 m . Owing to the SDSS observing strategy, stars 
brighter than go — 14 m will be saturated, and star counts 
will be not complete for magnitudes fainter than go = 22 m .2 
(Abazajian et al. I2004I I. Hence our work is limited for the 
magnitude range 15 m < go < 22 m for the evaluation of 
reliable Galactic model parameters. 

1 http://www.sdss.org/dr4/access/index.html 
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3 Density laws 

Disc structures are usually parametrized in cylindrical coor- 
dinates by radial and vertical exponentials, 

Di(x, z) = rii exp{-\z\/ Hi) exp(-(x - R )/hi), (2) 

where z = z Q +r sin(6), r is the distance to the object from 
the Sun, b is the Galactic latitude, z is the vertical distance 
of the Sun from the Galactic plane (24 pc, Juric et al. 2005 ), 
x is the planar distance from the Galactic centre, R$ is the 
solar distance from the Galactic centre (8 kpc, Reid 11993]) . 
Hi and hi are the scaleheight and scalelength, respectively, 
and rii is the normalized density at the solar radius. The 
suffix i takes the values 1 and 2 as long as the thin and thick 
discs are considered. 

The density law for the spheroid component is parame- 
terized in different forms. The most common is the de Vau- 
couleurs ( 1948 ) spheroid used to describe the surface bright- 
ness profile of elliptical galaxies. This profile has been de- 
projected into three dimensions by Young ( 1 19761 1 as 

D S (R) = n s exp[-7.669(i?/i? e ) 1 / 4 ]/(i?/i? e ) 7 / 8 , (3) 

where R is the (uncorrected) Galactocentric distance in spher- 
ical coordinates, R e is the effective radius and n s is the nor- 
malized local density. R has to be corrected for the axial 
ratio k = c/a, 

R=[x 2 + {z/n) 2 ] 1 / 2 , (4) 
where, 

z = r sin 6, (5) 

x= [Rl + r 2 cos 2 b - 2R Q rcosbcosl] 1/2 , (6) 
r being the distance along the line of sight and, b and I the 
Galactic latitude and longitude respectively, for the field un- 
der investigation. The form used by the Basel group is in- 
dependent of effective radius but is dependent on the solar 
distance from the Galactic centre: 

D S (R) = n s esp[10.093(l - i?/i? ) 1/4 ]/(i?/i? ) 7/8 - (7) 

Since star counts at high Galactic latitudes are not strongly 
related to the radial distribution, they are well suited to study 
the vertical distribution of the Galaxy. Hence, we adopted 
the scalelength of the thin and thick discs, 2.4 and 3.5 kpc 
respectively, of Juric et al. (2005) based on ~ 48 million 
stars, and we estimated the space densities for three com- 
ponents of the Galaxy at the solar distance, the scaleheights 
of the thin and thick discs, and the axial ratio of the halo 
by comparison of the space densities with the density laws 
in perpendicular direction. This is the procedure used by 
many researchers such as Phleps et al. (2000), Siegel et 
al. d2002b . Du et al. d2003l KBH (120041 Bilir, Karaali & 
Tu ncel ([2 005]), Phleps et al. J2005T) . BKG d2006b and Du et 
al. (127)061) . 



4 Metallicity calibration and absolute 
magnitude determination 

The metallicities of the sample stars with 0.12 < (g — r)o < 
0.95, which corresponds to 0.3 < (B — V)o < 1.1, were 
evaluated using the following calibration (KBT): 

[M/H] = 0.10 - 3.54^0.43 ~ 39.63^.43 + 63.51<^. 43 , (8) 

where <5o.43 is the normalized UV-excess in SDSS photom- 
etry. This parameter's range < 60.43 < 0.33 covers the 
metallicity interval -3 < [M/H] < 0.2 dex. 

Absolute magnitudes have been determined using the 
procedure described by KBT. They give the absolute mag- 
nitude offset from the Hyades main sequence as a function 
of both (3 — r)o colour and <5o.43 UV-excess, as follows: 

AMf = 03^.43 + C 2 S 2 43 + Ci(5o.43 + Cn, (9) 

where the coefficients Cj (i=0, 1, 2, 3) are functions of (g — 
r)o colour (Table 2) and are adopted from the work of KBT. 
AMg is defined as the difference in absolute magnitudes 
of a program star and a Hyades star of the same (3 — r)o 
colour: 

AM* = M* - Mf, (10) 

and the absolute magnitude for a Hyades star can be evalu- 
ated from the (Hyades) sequence, normalized by KBT: 

Mf = -2.0987( ff - rf - 0.0008(u - g) 2 
+0.0842( 5 - r)(u - g) + 7.7557(3 - r) 

-0.1556(u - 3) + 1.9714. (11) 

The KBT calibration is available for 4 < M(g) < 8 abso- 
lute magnitude interval where luminosity function is almost 
flat in the solar neighbourhood (Jahreiss & Wielen [T9971 ). 

In a conical magnitude-limited volume, the distance to 
which intrinsically bright stars are visible is larger than the 
distance to which intrinsically faint stars are visible. The 
effect of this is that brighter stars are statistically overrep- 
resented and the derived absolute magnitudes are too faint. 
This effect, known as Malmquist bias (Malmquist 1 1 9201 ), 
was formalized into the general formula, 

M(g) = M -a 2d -^l, (12) 

where M(g) is the assumed absolute magnitude, Mo is the 
absolute magnitude calculated for any star using KBT cali- 
bration, a is the dispersion of the KBT calibration, and A(g) 
is the differential counts evaluated at the apparent magni- 
tude 30 of any star. The dispersion in absolute magnitude 
calibration is around 0.25 mag, that produces a correction 
due to the Malmquist bias of less than 0.07 mag. The cor- 
rection of the Malmquist bias was applied to SDSS photo- 
metric data. Combination of the absolute magnitude M(g) 
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and the apparent magnitude go of a star gives its distance r 
relative to the Sun, i.e., 

[ 5 -M( ff )] = 51ogr-5. (13) 

Logarithmic space densities D* = log D + 10 have been 
calculated for the combination of three population compo- 
nents (thin and thick discs and halo), where D = N/AVi t 2\ 
AVi j2 = (7r/180) 2 (D/3)(r|-r?); □ denotes the size of the 
field (60 deg 2 ); r\ and r 2 denote the lower and upper limit- 
ing distances of the volume AV1.2; N is the number of stars 
per unit absolute magnitude; r* = [(rf + rf)^] 1 / 3 is the 
centroid distance of the volume A 14,25 and z* — r* sin(6), 
b being the Galactic latitude of the field center. The limiting 
distance of completeness is calculated from the following 
equation: 

[ 5 -M( 5 )] = 51ogn-5, (14) 

zi=rism(b) (15) 

where go is the limiting apparent magnitude (15 and 22, for 
the bright and faint stars, respectively), r; is the limiting 
distance of completeness, and M(g) is the corresponding 
absolute magnitude defining the interval (M\, M 2 ]. 

5 Estimation of the Galactic model 
parameters 

A x 2 method was employed to fit the analytical density laws 
given in Section 3 for combined all three population, i.e. 
thin and thick discs and halo, to the observational based 
space densities with the additional constraint of producing 
local densities consistent with those derived from Hipparcos 
(Jahreiss & Wielen 1997), a procedure applied in our pre- 
vious papers (KBH, BKG). All error estimates were made 
by changing Galactic model parameters until an increase or 
decrease by 1 in % 2 was achieved (Phleps et al. 2000). The 
comparison is shown in Fig. 3, and the resulting Galactic 
model parameters are given in Table 3 for two fields. Com- 
parison of the Galactic model parameters for two fields give 
interesting results: the scaleheight of the thin disc for two 
fields, 206(±9) pc (for north) and 198(±9) pc (for south), 
are rather close to each other. However this is not the case 
for the thick disc, i?=493(±12) and 579(±15) pc for the 
north and south fields, respectively. According to Narayan 

6 Jog (120021 1 and Momany et al. (2006) the difference be- 
tween the scaleheights is due to the warp and flare. The 
space density of the thin disc at the solar distance is almost 
the same for two fields, but there is a considerable differ- 
ence between the space densities of the thick discs, i.e. ~16 
and ~10 per cent relative the space density of the thin disc, 
for the north and south fields respectively. For the Galac- 
tic model parameters of the halo, one can say that there is 
a good agreement between the north and south ones. Espe- 
cially the axis ratio, n = 0.45 for both fields, is consistent 
with the values of Larsen ( 1 19961 k = 0.48), Chen et al. 
(20011 k = 0.55), and with the more recent one of Juric et 
al. (120051 k = 0.50). 



6 Metallicity 



We estimated the metal abundances of unevolved G type 
stars with 5 < M(g) < 6 by means of eq. (8), and we 
looked for a probable metallicity gradient in the vertical 
direction of the Galaxy. These are long lived stars which 
carry the heritage of the original chemical composition of 
the Galaxy, thus providing important information about the 
formation history of our Milky Way Galaxy. The metallic- 
ity distributions for 12 distance intervals are given in Fig. 
4 and Fig. 5 for the north and south fields, respectively. In 
the panels from (a) to (e) and from (i) to (1), in both fig- 
ures, one can notice a peak shifting towards low metallicies 
when one goes from short to large distances (r) relative to 
the Sun. Whereas there is a flat distribution for the panels 
(f), (g) and (h) which cover the distance interval 4 < r < 7 
kpc or 3.2 < z < 4.6 kpc, z being the distance from the 
Galactic plane. We adopted the abscissa of the peaks as the 
metal abundance of the corresponding distribution whereas 
for the flat distributions we used the medians. The metal 
abundances for 12 distance intervals mentioned above for 
north and south fields are given in Table 4. The number of 
stars, N, centroid distance, r* , and the corresponding dis- 
tance in the vertical direction, z*, are also given in the same 
table. 

The metallicity versus z* diagrams for two fields are 
given in Fig. 6. Three substructures are conspicuous in this 
figure, i.e. the smooth tendencies for the intervals z* < 3 
kpc and 5 < z* < 10 kpc, and the steep one for the in- 
terval 3 < z* < 5 kpc. The tendencies in the first inter- 
val correspond to low metallicity gradients, d[M/H]/dz = 
-0.16(±0.02), and -0.22(±0.02) dex kpc" 1 for the north 
and south fields, respectively. The interval z* < 3 kpc in- 
volves the thin disc which is formed by dissipative collapse 
of the corresponding part of the Galaxy. Hence, the cited 
metallicity gradients are (partly) due to the formation proce- 
dure of the thin disc. The metallicity gradients detected for 
two fields for the interval 5 < z* < 10 kpc are almost zero, 
d[M/H]/dz = -0.036(±0.001) and -0.034(±0.001) dex 
kpc -1 for the north and south fields, again consistent with 
the formation scenarios of the halo (see Section 7 for de- 
tail). The metallicity gradients for two fields for the inter- 
val 3 < z* < 5 kpc are rather different from each other, 
i.e. d[M/H]/dz = -0.36(±0.12) and -0.18(±0.01) dex 
kpc -1 , for the north and south fields respectively. Also, the 
first value is the highest gradient detected in this work. The 
interval 3 < z* < 5 kpc covers both thick disc and halo 
stars. Hence, the metallicity gradients cited for this inter- 
val probably correspond to the transition from thick disc to 
halo. The combination of the data in Table 4, for two fields 
individually, give metallicity gradients almost equal to each 
other for the whole Galaxy, d[M/H]/dz = -0.29(±0.04) 
and — 0.30(±0.04) dex kpc -1 for the north and south fields 
respectively. 
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7 Conclusion 

We estimated Galactic model parameters for two fields in 
the anti-centre direction of the Galaxy (7 = 180°) and sym- 
metric relative to the Galactic plane, b = +45° (north field) 
and b = —45° (south field). The size of the fields are large, 
avoiding any fluctuation in the values of the model param- 
eters. We adopted the scalelengths for thin and thick discs 
from Juric et al. (120051 1. h\ = 2.4 and h% = 3.5 kpc respec- 
tively, determined by means of ~48 million stars, and we 
estimated the space densities at the solar distance for three 
Galactic components, scaleheights of thin and thick discs, 
and the axis ratio for the halo. 

No conspicuous differences could be detected between 
the corresponding Galactic model parameters for the thin 
discs of two fields. However this is not the case for the thick 
discs. Actually, the scaleheights are 493 and 579 pc, and the 
space densities at the solar distance normalized to the thin 
disc are 16 and 10 per cent for the north and south fields 
respectively. According to Narayan & Jog (2002), and Mo- 
many et al. (120061 1. these differences originate from warps 
and flares. For halo, the axis ratio of two fields are exactly 
the same, k = 0.45, and the space densities at the solar dis- 
tance are rather close to each other. 

Comparison of our model parameters with the ones of 
Juric et al. ( 120051 1 from which the scalelengths of the thin 
and thick discs were adopted shows both agreement and dif- 
ferences between two sets of data. For example, Juric et al. 
(2005) give n ~ 0.5 for the axis ratio of the halo which is 
almost equal to the value cited in our work. However, their 
scaleheights of the thin and thick discs are different than 
ours, especially the scaleheights of the thick disc for both 
fields are less than the half of their scaleheight, i.e. H2 ~ 
1200 pc. The scaleheight of the thick disc cited by Juric et 
al. d2005) reminded us the original value of Gilmore & Reid 
(1983), H =1-1.4 kpc, however it is rather different than 
the updated value (cf. Chen et al. 2001 ). Another conspicu- 
ous difference is between the normalized space densities of 
the thick disc at the solar distance. Juric et al. (2005 ) give 
a local thick-to-thin disc normalization 4 per cent which is 
0.25 and 0.40 times of the values cited in our work for the 
north and south fields respectively. Our local densities for 
the thick disc are close to the updated ones (cf. Chen et al. 
120011) whereas the one of Juric et al. (2005 ) is close to the 
original one, ~ 2 per cent. Despite these differences, the 
Galactic model parameters cited in our work are consistent 
with the corresponding ones appeared in the literature (see 
Table 1) with two exceptions, i.e. the local space density 
normalization for the thick disc for the north field, 16 per 
cent, and the scaleheight of the thick disc for the north field, 
493 pc. 

Three substructures can be revealed for the metallicity 
distribution of stars with 5 < M(g) < 6: 
(1) For z* < 3 kpc the metallicity gradients are smooth 

and close to each other for two fields, d[M/H]/dz = 

-0.16(±0.02), and -0.22(±0.02) dex kpc" 1 for the 

north and south fields respectively. 



(2) For the interval with stars at large distances, i.e. 5 < 
z* < 10 kpc, the metallicity gradient is d[M/H]/dz ~ 
-0.03(±0.004) dex kpc" 1 for two fields. 

(3) For the intermediate distance interval 3 < z* < 5 kpc, 
the metallicity gradients for two fields are rather dif- 
ferent from each other, d[M/H]/dz = -0.36(±0.12) 
and — 0.18(±0.01) dex kpc -1 for the north and south 
fields respectively. Also, the tendency of the metallicity 
for stars in the north field is rather steep. The combi- 
nation of the metallicity distributions in these substruc- 
tures gives a metallicity gradient of d[M/H]/dz ~ —0.30 
(±0.04) dex kpc -1 for two fields. 

The result cited in (1) is consistent with the formation 
scenario of the thin disc, i.e. this component of the Galaxy 
was formed by dissipative collapse. The smooth gradient 
is due to the contamination of the thick disc for which no 
metallicity gradient is expected. The interval 5 < z* < 10 
kpc covers predominantly the halo stars, and the metallicity 
close to zero is the expected one. That is, halo is subject to 
some mergers or accreation of the objects formerly formed 
outside of our Galaxy. The most interesting result is for the 
interval 3 < z* < 5 kpc: the metallicity gradients for the 
north and south fields are different from each other. This 
interval is partly occupied by the thick disc stars. Hence, the 
metallicity gradient cited for this interval corresponds to the 
transition region from the thick disc to the halo. We remind 
the reader the differences between the scaleheights and local 
space densities of the north and south fields cited above. 
These findings encariage us to argue that the structure of 
the thick disc is Galactic latitude dependent. 

The metallicity gradients cited in our work are consis- 
tent with the corresponding ones appeared in the literature. 
For example, Trefzger et al. ( 1995) estimated a metallicity 
gradient of d[M/H]/dz = - 0.23(± 0.04) dex kpc" 1 for 
z < 4 kpc, and Karaali et al. ( 120031) give d[M/H]/dz = 
-0.20, -0.10, and 0.20 dex kpc -1 for the intervals z < 
5, 5 < z < 8, and z < 8 kpc. For a more recent cita- 
tion, we can give that of Du et al. (120041 ): d[M/H]/dz = 
-0.37(±0.01), -0.06(±0.01) and -0.17(±0.04) dexkpc" 1 
for z < 4, z > 4 and z < 15 kpc. 
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Fig. 3 Comparison of observed (symbols) and calculated (lines) space density functions combined for stars of all three 
population components and with absolute magnitudes 4 < M(g) < 8 for two star fields. 



© 2006 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



Astron. Nachr. / AN (2006) 



9 



0) 

cu 

-I— ' 

w 



_g 100 




Fig. 4 Metallicity distributions for stars with 5 < M(g) < 6 absolute magnitude within the r-distance ranges (in kpc) 
for north field: (a) (0.5,1.5], (b) (1.5,2], (c) (2,2.5], (d) (2.5,3], (e) (3,4], (f) (4,5], (g) (5,6], (h) (6,7], (i) (7,8], (j) (8,9], 
(k) (9,12] and (1) (12,15]. The curves are Gaussian fits with mean metallicities and dispersion in Table 3. The arrow in the 
panels shows the median. 



Table 1 Previous Galactic models. Symbols: TN denotes the thin disc, TK denotes the thick disc, S denotes the spheroid 
(halo), R e is the effective radius and n = c/a is the axis ratio. The figures in the parentheses for Siegel et al. (120021) are the 
corrected values for binarism. The asterisk denotes the power-law index replacing R e . 



H (TN) 


h (TN) 


n (TK) 


H (TK) 


h (TK) 


n(S) 


Re(S) 


ft 


Reference 


(pc) 


(kpc) 




(kpc) 


(kpc) 




(kpc) 






310-325 




0.0125-0.025 


1.92-2.39 










Yoshii d 1982b 


300 




0.02 


1.45 




0.0020 


3.0 


0.85 


Gilmore & Reid 11983b 


325 




0.02 


1.3 




0.0020 


3.0 


0.85 


Gilmore 11984b 


280 




0.0028 


1.9 




0.0012 






Tritton & Morton 119841 


200-475 




0.016 


1.18-2.21 




0.0016 




0.80 


Robin &Creze 11986b 


300 




0.02 


1.0 




0.0010 




0.85 


del Rio &Fenkart 11987b 


285 




0.015 


1.3-1.5 




0.0020 


2.36 


flat 


Fenkartetal. 11987b 


325 




0.0224 


0.95 




0.0010 


2.9 


0.90 


Yoshii et al. 11987b 


249 




0.041 


1.0 




0.0020 


3.0 


0.85 


Kuijken & Gilmore |l989|l 


350 


3.8 


0.019 


0.9 


3.8 


0.0011 


2.7 


0.84 


Yamagata & Yoshii 11992b 


290 






0.86 






4.0 




von Hippel & Bothun |Tl993l 


325 




0.0225 


1.5 




0.0015 


3.5 


0.80 


Reid &Majewski 11993b 


325 


3.2 


0.019 


0.98 


4.3 


0.0024 


3.3 


0.48 


Larsen i 1996b 


250-270 


2.5 


0.056 


0.76 


2.8 


0.0015 


2.44-2.75* 


0.60-0.85 


Robin et al. 11996b; Robin et al. 12000b 


290 


4.0 


0.059 


0.91 


3.0 


0.0005 


2.69 


0.84 


Buser, Rong & Karaali 11998111999b 


240 


2.5 


0.061 


0.79 


2.8 






0.60-0.85 


Ojhaetal. 11999b 


330 


2.25 


0.065-0.13 


0.58-0.75 


3.5 


0.0013 




0.55 


Chenetal. 1200 lb 


280 


2-2.5 


0.06-0.10 


0.7-1.0 


3-4 


0.0015 




0.50-0.70 


Siegel et al. 12002b 


(350) 


2-2.5 


0.06-0.10 


(0.9-1.2) 


3-4 


0.0015 




0.50-0.70 


Siegel et al. d2002t 


320 




0.07 


0.64 




0.0013 




0.58 


Du et al. 120031 


280 


2.4 


0.04 


1.2 


3.5 






0.50 


Juric et al. J20051 
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Fig. 5 Metallicity distributions for stars with 5 < M(g) < 6 absolute magnitude within the r-distance ranges (in kpc) 
for south field: Symbols are the same in Fig. 4. 



Table 2 Numerical values for the coefficients a (i=0, 1, 2, 3) in eq. (9), adopted from the work of KBT (2005). 



(.9 - r )o 


C3 


C2 


Cl 


CO 


(0.12,0.22] 


-68.1210 


26.2746 


2.2277 


-0.0177 


(0.22,0.32] 


-32.5618 


6.1310 


5.7587 


0.0022 


(0.32,0.43] 


8.2789 


-7.9259 


6.9140 


0.0134 


(0.43,0.53] 


-23.6455 


-0.4971 


6.4561 


0.0153 


(0.53,0.64] 


0.2221 


-5.9610 


5.9316 


-0.0144 


(0.64,0.74] 


-47.7038 


0.1828 


4.4258 


-0.0203 


(0.74,0.85] 


-52.8605 


12.0213 


2.6025 


0.0051 


(0.85,0.95] 


-15.6712 


7.0498 


1.6227 


-0.0047 



Table 3 Galactic model parameter values estimated for star fields. The columns give: logarithmic local space density n*\ 
scaleheight H; local space densities relative to thin disc, for thick disc nijnx and for halo ns/ni; axis ratio of halo k. 



Field 


Parameter 


Thin disc 


Thick disc 


Halo 


North 


n* 


7.40 ± 0.09 


6.61 ±0.03 


4.90 ±0.06 




H(pc) 


206 ±9 


493 ± 12 






n/ni (per cent) 




16.03 ± 3 


0.31 ± 0.03 




K 






0.45 ± 0.06 


South 


n* 


7.43 ± 0.09 


6.44 ± 0.03 


4.79 ± 0.09 




H(pc) 


198 ±9 


579 ± 15 






n/ni (per cent) 




10.40 ± 3 


0.23 ± 0.03 




K 






0.45 ± 0.09 
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Fig. 6 Mean metal abundances versus z*-distance for 12 z* intervals, showing the metallicity gradients — 0.29(±0.04) 
and -0.30(±0.04)dexkpc- 1 up to 10 kpc for overall distribution, -0.16(±0.02) and -0.22(±0.02) dexkpc" 1 forz* < 3 
kpc, -0.36(±0.12) and -0.18(±0.01) dex kpc" 1 for 3 < z* < 5 kpc, -0.036(±0.001) and -0.034(±0.001) dex kpc" 1 
for 5 < z* < 10 kpc for north and south fields, respectively. 
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Table 4 Metallicity distribution for 12 distance intervals, r* being the centroid distance, z* centroid distance to the 
Galactic plane, [M/ H] mean metallicity and N number of stars. 









North 




South 




n - r 2 


r* 


2* 


[M/H] 


N 


[M/H] 


N 


(kpc) 


(kpc) 


(kpc) 


(dex) 




(dex) 




0.5-1.5 


1.20 


0.85 


-0.39 ± 0.05 


2988 


-0.34 ±0.05 


3109 


1.5-2.0 


1.78 


1.26 


-0.41 ±0.07 


2080 


-0.45 ±0.09 


2129 


2.0-2.5 


2.28 


1.61 


-0.47 ±0.05 


1666 


-0.56 ±0.12 


1684 


2.5-3.0 


2.77 


1.96 


-0.52 ±0.04 


1313 


-0.63 ±0.09 


1291 


3.0-4.0 


3.57 


2.52 


-0.66 ±0.08 


1904 


-0.71 ±0.16 


1678 


4.0-5.0 


4.56 


3.22 


-1.28 ±0.90 


1323 


-1.31 ±0.91 


975 


5.0-6.0 


5.54 


3.92 


-1.68 ±0.87 


1026 


-1.44 ±0.92 


635 


6.0-7.0 


6.54 


4.62 


-1.79 ±0.87 


801 


-2.56 ±0.90 


460 


7.0-8.0 


7.53 


5.33 


-2.26 ±0.21 


700 


-2.40 ±0.34 


375 


8.0-9.0 


8.53 


6.03 


-2.29 ±0.31 


581 


-2.43 ±0.33 


360 


9.0-12.0 


10.71 


7.57 


-2.34 ±0.55 


1167 


-2.48 ±0.40 


853 


12.0-15.0 


13.66 


9.66 


-2.42 ±0.40 


716 


-2.55 ±0.56 


660 
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